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ABSTRACT 



We report on Chandra observations of the black widow pulsar, PSR B1957+20. Evidence 
for a binary-phase dependence of the X-ray emission from the pulsar is found with a deep ob- 
servation. The binary-phase resolved spectral analysis reveals non-thermal X-ray emission of 
PSR B1957+20, confirming the results of previous studies. This suggests that the X-rays are 
mostly due to intra-binary shock emission which is strongest when the pulsar wind interacts with 
the ablated material from the companion star. The geometry of the peak emission is determined 
in our study. The marginal softening of the spectrum of the non-thermal X-ray tail may indicate 
that particles injected at the termination shock is dominated by synchrotron cooling. 

Subject headings: binaries: eclipsing — stars: individual (PSR B1957+20) — stars: neutron 



1. Introduction 

The widely accepted scenario for the for- 
mation of a millisecond pulsar (MSP) is that 
an old neutron star has been spun up to mil- 
lisecond periods in a past accretion phase by 
mass and angular momentu m transfer f r om a 
binary late-type companion (Al par et al. 1982L 
Radhakri shnan &: Srinivasa Once the 

accretion has stopped, the relativistic magne- 
tized pulsar wind, which is believed to carry 
away the pulsar rotational energy and angular 
momentum, may be able to ablate and eventu - 
ally evaporate its companion ( Rasio et al. 19891) . 
Close binary systems with MSPs are a sub- 
ject of special interest since they are thought 
to be the missing link between low-mass X- 
ray binaries (LMXBs) and isolated MSPs. Dis- 
coveries of the eclipsing binary pulsar systems , 
such as PSR B19 57+20. dFruchter et al.l Il988bh. 
PSR J2051-0 827 (IStappers et al.lll996l) a nd PSR 
J1023+0038 (jArchibald et al.l 120091 l2010h in the 
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Galactic plane and 47 Tuc W (PSR J0024 



7204W: 



Bogdanov et"afl 120051) and PSR J1740- 



5340 (jD'Amico et al.l 120011 ) in globular clusters, 
gave support to this formation scenario. Studying 
these binary systems provides a wealth of infor- 
mation not only on the evolutionary history of 
isolated MSPs but also on the physical details of 
the pulsar's high energy emission properties. 

PS R B1957+20 was disc overed at Arecibo in 
1988 (|Fruchter et al.lll988bl) . It is in a binary 



system with a 0.025 M Q companion in a 9.16-hr 
orbital period. The pulsar has a spin period of 
1.6 ms, the third shortest among all known MSPs. 
Its period derivative of P = 7.85 x 10~ 21 s s -1 
implies a spin-down energy of E = 7.48 ± 3.61 x 
10 34 erg s _1 , a characteristic spin-down age of 
> 3.4 x 10 9 yrs, and a d ipole surface magneti c 
field ofS ± = 1.12 x 10 8 G (|Guillemot et alj|2012ft . 
For a radio dis persion measure inferr ed distance 
of 2.5 ± 1.0 kpc (|Cordes fc Lazioll2002f ). the pulsar 
moves through the sky with a supersoni c velocity 
of 220 km s -1 (jArzoumanian et al.lll994[ ). The in- 
teraction of a relativistic wind flowing away from 
the pulsar with the interstellar medium (ISM) 
produces an Ha bow shock nebula which was the 
first nebula to be found arou nd a "recycled" pul- 
sar (jKulkarni fc Hesterlll988l ). For approximately 
10% of this orbit, the radio emission at a frequency 
of 430 MHz from the pulsar is eclipsed by mate- 
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rial ablated by the 



ion dFruchtcr et alJUMSSO) . U ptical observation; 
by Fruchter et al.1 ( 1988at ) and Ivan Paradiis et al 



pulsar wind from the compan 



19883). Optical observations 



(|1988I ) revealed that the pulsar wind consisting of 
electromagnetic radiation and high-energy parti- 
cles is ablating and evaporating its white dwarf 
companion star. This rarely observed property 
gave the pulsar the name black widow pulsar. Sub- 
sequent analyses o f the optical ligh t curv e for the 
binary system by iRevnolds et al. ( 2007 ) gave a 
constraint to the system inclination of 65° ± 2° 
for a pulsar in the mass range 1.3 — 1.9 Mq. The 
effective temperatures of T — 2900 ± 110 K and 
T = 8300 ± 200 K for the unilluminatcd side and 
the illuminated side of the companion were also 
obtained in their studies. 

Several X-ray studies of PSR B195 7+20 have 
been carried out in the past two decades ( Kulkarni et al 
19921 : IStappers et albOoilHuang fc Beckerl2007l ). 
The X-ray emission of PSR B1957+20 is found 
to be non-thermal dominated and best modeled 
with a single power-law spectrum, which indi- 
cates that the X-rays originate from the shock 
interaction of the pulsar wind with the wind of 
the comp anion star or from t h e pulsar magne- 
tosph ere (jStappers et al.l l2003t iHuang fc Becker 
20071 ). Searching for an X-ray pulsation at the 
radio pulsar's rotation period and a modulation 
of the X-ray emission from the PSR B1957+20 
system with orbital phase is another important 
subject, which may help to discriminate the ori- 
gin of the X-ray emission. Recently a 4 — a 
detection of X-ray cohere n t pul sation was re- 
ported bv iGuillemot et al.1 (|2012l ). In addition, 
Huang fc Becked (|2007h found a strong correlation 
of the pulsar's X-ray flux with its orbital period. 
However, due to the short exposure we could not 
know whether the flux modulation was periodic 
and given the limited photon statistics it was not 
possible to investigate any spectral variation as a 
function of orbit phase or to determine the exact 
geometry of the peak emission. Repeated coverage 
of the binary orbit in a longer Chandra observa- 
tion would provide us a better photon statistic 
and allow us to determine the emission geometry 
with higher accuracy. 

In this paper, we report on an archival Chandra 
observation of the black widow pulsar. This inves- 
tigation offers further insight into the properties of 



this remarkable binary. In §2, we summarize the 
details of the Chandra observation and the data 
reduction. A search for the X-ray orbital modula- 
tion from this system was described in §3. In §4 
and §5, we present the results of the X-ray spa- 
tial and spectral analysis. Finally, we discuss the 
physical implications of the observed results in §6. 

2. Observation 

A Chandra observation aimed on PSR B1957+20 
was performed on 2008 August 15 (ObsID 9088) 
using the back-illuminated chip ACIS-S3 with an 
uninterrupted 169-ks exposure. The data was 
configured in the VFAINT telemetry mode. Data 
reduction and analysis were processed with Chan- 
dra Interactive Analysis Observations (CIAO) ver- 
sion 4.3 software and the Chandra Calibration 
Database (CALDB) version 4.4.1. The level 1 
data with background cleaning were used in our 
study. Data analysis was restricted to the energy 
range of 0.3 - 8.0 kcV. 

For the timing and spectral analyses of the 
black widow pulsar, we extracted the photons from 
a circular region centered at the radio timing posi- 
tiorQ, RA(J2000)=19' , 59 m 36 1 r77, Dec=20°48'15'/12, 
with a radius of 2" which encloses 90% of the total 
source energy at 1.5 keV. 

3. Timing Analysis 

For the timing analysis, we first extracted the 
photons from the aforementioned circle and trans- 
lated the photon arrival times to the solar system 
barycenter by using the CIAO tool axbary. The 
JPL DE200 solar system ephemeris was used for 
the barycentric correction to ensure consistency 
with the radio ephemeris. We note that search- 
ing for the X-ray pulses at the spin period of 
PSR B1957+20 was precluded by the inappropri- 
ate temporal resolution of this Chandra observa- 
tion with a frame time of 3.2 s. 

As the Chandra observation covers over five 
consecutive binary orbits, by plotting a light curve 
of the X-ray source counts versus the orbital phase 
(see Figure Q] left panel) we can confirm that the 
X-ray flux is not steady with time. We also applied 
a Kolmogorov-Smirnov (KS) test to the unbinned 



ifrom the ATNF Pulsar Catalogue 
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light curve data in order to have a bin-independent 
statistical evaluation of the X-ray emission vari- 
ability. Calculating the corresponding KS proba- 
bilities between our data set and the cumulative 
distribution function generated by assuming a con- 
stant X-ray flux indicated a marginal deviation 
between these two distributions. The significance 
for an intra-orbital flux modulation from this test 
is only at the level of ~81%. A more significant 
result, i.e. ^96%, is obtained if we restrict the 
analysis to the soft energy band below 2 keV. 

In order to search for a modulation of the X-ray 
flux as a function of orbital phase, we first selected 
X-ray data covering 5 complete and consecutive 
orbits and then used the radio timing ephemeris 
of PSR B1957+20 from a pulsar catalog provided 
by Lucas Guillemot to fold a light curve at the 
orbital period (see Figure Q] right panel). Using a 
X 2 -test, the significance for a flux modulation over 
the observed orbit was found to be ~ 99%. 

4. Spatial Analysis 

Figure [2] (left panel) shows the Chandra ACIS- 
S3 image in the energy band 0.3-8 keV of the 
field around PSR B1957+20. This image was cre- 
ated by using an adaptive smoothing algorithm 
with a Gaussian kernel of a < 3 pixels in or- 
der to probe the detailed structure of faint diffuse 
emission. Both the pulsar and an extended X-ray 
feature (hereafter the "tail" ) , protruding from the 
pulsar position, can be clearly seen in this image. 
The length of the tail with its orientation to the 
northeast is about 25 arcsec. The right panel in 
Figure [2] presents the Ha image overlaid with the 
X-ray contours. This Ha image is obtained from 
Taurus Tunable Filter service mode observations 
on the Anglo Australian Telescope in 2000. An Ha 
bow-shock emission aligned with the proper mo- 
tion direction of the PSR B1957+20 system fades 
to the background in ~ 40 ". Comparing with th e 



previous result reported bv lStappers et al.l (|2003h . 



we found that the X-ray tail is more extended due 
to the longer exposure time of this Chandra ob- 
servation. 



2 ftp: / / www.cenbg.in2p3.fr / astropart /lucas / report /1959+2048.html 



5. Spectral Analysis 
5.1. PSR B1957+20 

We used the CIAO tool dmextract to extract 
spectra of the source and two nearby source-free 
background regions. Response files were con- 
structed by using the CIAO tool mkacisrmf and 
mkarf. The extracted spectra were binned with 
at least 30 source counts per bin. Background- 
subtracted spectral modeling was performed with 
XSPEC (version: 12.7.0) using data in the energy 
band 0.3-8.0 keV. 

Assuming that the X-ray emission originates 
from the in tra-binary shock or the pulsar's mag- 
netosphere (IStappers et al.ll2003 ; Huang fc Becker 
120071 : iGuillemot et al.l 2012 ). we expect the radi- 
ation to be synchrotron. To test this hypothesis, 
we fitted the spectrum with an absorbed power- 
law model (PL). Unexpectedly, a single PL model 
cannot provide any statistically acceptable de- 
scription of the observed spectrum (i.e. xt > 2). 
We then tested whether a single blackbody (BB) 
model, a double PL model, or a composite spectral 
model consisting of a PL and a BB can provide an 
appropriate modeling of the data. Neither the one- 
component BB model nor the double-component 
PL model yields any physically acceptable fits. In- 
stead, the composite model gives a better descrip- 
tion of the X-ray spectrum. The inferred tempera- 
ture of 0.18 keV and 0.30 keV is much too high for 
the cooling surface of the neutron star, but may 
be explained by the heated po lar cap region (c.f. 
Zavlinll2007l:lTakata et aTll2012l ). Considering that 



a portion of the X-ray emission could originate 
from a thermal plasma within or around the bi- 
nary system, possibly from the active corona of the 
secondary star or the plasma respon sible for the 
radio eclipses ( Bogdanov et al.ll2011 ). a MEKAL, 
a thermal bremsstrahlung (TB), a MEKAL+PL, 
and a PL+TB model were tested in our study. We 
found those models cannot provide a better or an 
acceptable description of the X-ray spectrum of 
the pulsar system. Besides, we also examined the 
X-ray spectrum with more complicated spectral 
models, e.g. PL+PL+BB or PL+BB+BB. Al- 
though an acceptable fit can be yielded by a three- 
component model, i.e. PL+BB+BB, no physical 
interpretation can be appropriately applied to this 
binary system. Therefore, we suspect a depen- 
dence of the X-ray spectrum of PSR B1957+20 on 
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Fig. 1. — Left: X-ray emission from PSR B1957+20 in the energy range of 0.3-8.0 keV as a function of the 
pulsars orbital phase ((f). Right: A folded light curve at the orbital period. Two orbital cycles are shown for 
clarity. The background noise level is found to be at ~ 0.1 counts/bin. The phase zero (4>= 0.0) corresponds 
to the ascending node of the pulsar orbit. Error bars indicate the la uncertainty. The blue shaded regions 
between the orbital phases 0.2-0.3 and 1.2-1.3 mark the radio eclipse of the black widow pulsar. The 
phase-resolved spectrum covering the eclipsing region was extracted from the gray shaded regions (see §5.1). 
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Fig. 2. — Left panel: Chandra ACIS-S3 image in the energy band 0.3-8 keV of the black widow pulsar system 
smoothed with an adaptive Gaussian filter. The green circle with the 2.0" radius indicates the source region 
we used in this study while two segments of the X-ray tail are chosen from the red rectangular regions. Right 
panel: The Ha image taken from the Anglo Australian Telescope is overlaid with the X-ray contours. The 
green contour levels are shown at 0.4, 0.8, 3.0, 15.4, and 84.8% of the peak x-ray surface brightness. The red 
cross indicates the radio timing position of PSR B1957+20. The optical residuals correspond to incompletely 
subtracted stars. 
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its orbital phase due to the variability observed in 
its X-ray flux level which seems to correlate with 
its orbital period. 

To investigate whether the X-ray spectral be- 
havior of PSR B1957+20 varies across the orbit, 
we analysed the X-ray spectra within the orbital 
phase of (f> = 0.05 — 0.45 which covers the ingress, 
eclipsing, and egress region and outside the afore- 
mentioned region (cf> = 0.45 — 1.05) separately. 
The source and background spectra were extracted 
from the same circles as mentioned before and dy- 
namically binned in accordance with the photon 
statistic in each dataset. The net count rates for 
the spectral analysis inside and outside the eclips- 
ing region arc 7.49 ± 0.33 and 5.88 ± 0.25 counts 
per kilosecond, respectively. We found that the 
binary-phase resolved spectral analysis reveals a 
non-thermal emission nature of the detected X- 
rays and each of the observed spectra can be well 
described by a single PL model with different 
photon indices, which indicates that its spectral 
behavior is orbital dependent. 

The spectral parameters inferred from the PL 
fits as well as the best-fit hydrogen column density 
Nh are found to be c onsist ent wi th those reported 



(120071) 



by IStappers et alj (|2003|) and iHuang fe Becker 
For comparison, a fixed Nh value of 



1.8 x 10 21 enr 2 obtained bv lStappers et all (|2003l ) 
is also applied in our spectral analysis. The results 
of the spectral fits are summarized in Table 1. 

5.2. The X-ray Tail 



Stappers et al.l (|2003[ ) uncovered an X-ray tail 
with a position angle opposite to the proper mo- 
tion direction of PSR B1957+20 with a Chandra 
observation taken in 2001. However, with low 
photon statistics, a detailed spectral analysis was 
restricted. The much longer exposure time of the 
latest archival Chandra observations allows us to 
study the X-ray tail in greater detail. For com- 
parison, we first selected a box of 16" x 6" with 
an orientation along the proper motion direction 
as the region of the X-ray tail of PSR B1957+20. 
The background-subtracted count of 108 ± 13 was 
collected for the spectral analysis. With higher 
photon statistics, we found an absorbed single PL 
model fits the X-ray spectrum of the tail well, 
which implies that the X-ray emission originates 
from the pulsars interaction with the ambient 



medium. Its non-thermal X-ray spectrum points 
to synchrotron emission from energetic particles 
from the pulsar wind. 

A softening of the spectrum of the X-ray tail 
as a function of the distance from the pulsar is 
expected if synchrotron cooling of the particles in- 
jected at the termination shock is dominated. For 
the purpose of investigating the possible spectral 
variation, we performed a spatially-resolved spec- 
tral analysis using two separate extraction regions 
along the tail. We refer the segment close to the 
pulsar as region 1 and the further one as region 2 in 
the following. The background-subtracted counts 
are 68 ± 14 for region 1 and 29 ± 12 for region 2, 
respectively. In order to better constrain the spec- 
tra l properties, we f o llowed the method adopted 
by IJohnson fc Wand (|2010f ) and iHui et al.l (|2012l ) 
to jointly fit individual power-law models for the 
X-ray spectra of these two regions assuming the 
column density, Nh, does not change significantly 
along the proper motion direction. Due to low 
photon statistics in each segment, we then fixed 



the N h value at 1.8 x 10 cm 2 (jStappers et al 



20031 ) in the joint fit. The best- fit spectral param- 
eters are shown in Table 1. An indication for such 
a spectral variation was found in this study. 

6. Summary & Conclusion 

We have searched for the orbital modula- 
tion of the X-ray emission from PSR B1957+20. 
Analysing this data set with a x 2_ test and a 
Kolmogorov-Smirnov test revealed a marginal 
intra-orbital flux modulation, which suggests that 
the non-thermal X-rays from PSR B1957+20 are 
mostly due to intra-shock emission at the inter- 
face between the pulsar wind and the ablated 
material from the companion star. The pulsar 
wind electrons and positrons are accelerated and 
randomized by the shock and emit the X-rays 
via the synchrotron process. Such shock emis- 
sion has also been suggested to explain the vari- 
able X-ra y flux as a function of orbital phase in 
47 Tuc W dBogdanov et al.ll2005h and P SR, J1740- 
5340 in NGC 6397 (IHuang fc Beckeil l2010f) and 
PSR J1023+0038 (lArchibald et all I2009L hoiot 
Tarn et al.ll201oHBogdanov et al.ll201ll) . 

In Figure [IJ the observed flux peaks just 
before and after the pulsar eclipse can be in- 
terpreted as the Doppler effect caused by the 
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Table 1 



X-RAY SPECTRAL PARAMETERS OF THE PSR B1957+ 20 BINARY SYSTEM 



Orbital Phase 


Model 


N H 

(10 21 cm" 2 ) 




r/kT 

/(keV) 




F x a 

(10 -14 ergs cm -2 s^ 1 ) 


Xi/d.ol 


Pulsar 


0.0-1.0 


PL 

PL 
PL+BB 
PL+BB 


0.83t H 
1.80 b 

< 0.61 

1.80 b 


r= 
r= 


1.96 + 0.12 
2.26 c 

lmtoil, kT=0.30; 
lmtlil, kT=o.i8; 


f0.02 
-0.05 
f0.02 
-0.03 


6 D5+ - 26 

D.UO_ 2 8 

7 99+O.I3 
' - zz -0.14 

5.39l°™ 

7 i O+0.12 


1.86/18 
2.17/19 
1.48/16 
1.75/17 


0.05-0.45 


PL 
PL 


1.80 b 




1.80±° ;1 9 7 
2.09 + 0.10 




D nc+0.53 
°- uo -0.70 
90+O.43 
J - ZO -0.39 


0.75/13 
0.92/14 


0.45-0.05 


PL 
PL 


89+ - 54 

u - oa -0.50 

1.80 b 




9 -, 7 +0.21 

z - x ' -0.20 

2 49+ 11 
z -^ y -o.io 




6 12+ 036 

u - lz -0.38 

7 61 +0 ' 24 

' -OJ-_0.19 


0.99/15 
1.10/16 



Segment 


Model 


N H 

(10 21 cm- 2 ) 


r 


F x a 

(1Q- 14 ergs cm~ 2 s" 1 ) 


X^/d.o.f 


Tail 


Whole tail 


PL 
PL 


9 77+3.21 

Z.I /_ 2 .42 

1.80 b 


9 fi 9 + 1.23 
z - DZ -0.84 

2 30+ 39 

Z ' OU -0.36 


0.891°;}° 
0.7lig;S? 


0.77/4 
0.64/5 


Region 1 
Region 2 


PL 
PL 


1.80 b 


1 t-7+0.31 
i -°'-0.26 

2 14+ - 46 

Z ' 1 ^-0.37 


46+ 013 
u - w -o.io 

S3+ 016 


1.28/18 





a Unabsorbed X-ray flux in the energy range of 0.3-8.0 keV. 
b Thc hydrogen column density Nh is fixed at 1.8 x 10 21 cm -2 . 

c No error was calculated since the reduced \ 2 is larger than the maximum value of 2. 
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bulk flow in the down stream region. If the 
post-shocked wind flows toward (or away from) 
the Earth, the Dopplcr effect increases (or de- 
creases) the observed flux from the flux for an 
isotropic case. The shock geometry is controlled 
by the ratio of the momentum fluxe s of the pul 



sar wind to the stellar w ind (c.f. ICanto et al 
119961 : lAntokhin et al.lliool . For PSR B1957+20, 
the ob served orbital perio d derivative Pb ~ 
1CT 11 (|Fruchter et all Il990h suggests the mass 
loss rate M* ~ M^P^ 1 ~ 10" 10 M Q yr" 1 , 
where M* = 0.02 M Q and P b = 33001 s is 
the orbital period. With the mass-loss rate of 
M ~ 10" 10 M© yr" 1 , the ratio of the momen- 
tum fluxes of the pulsar wind and the stellar 
wind is in the order of 77 = L s d/ M*v es c ~ 



10(W10 35 erg S - 1 )(M*/1O- 1O M yr^X^/SOO km s" 
indicating the companion star is confined by the 
pulsar wind and the shock. With momentum ratio 
X) ~ 10, the ope ning angle of the cone-like shock is 
~ 50 - 60° (c.f. lEichler fc Usovll 19931 ICanto et al ' 



where the observed fluxes arc maximum and min- 
imum, respectively. For example, the observed 
ratio F max /F m i n ^2 — 3 and the photon index 
a ~ 2 implies the flow velocity /3/ ~ 0.12 — 0.19 for 
Of, max = 4 5° and fl, ~ 0.2- 0.32 for f , max = 65°. 
Note that iGuillemot et all (|2012l) estimates the 
fraction of the pulsed emissions as about 30 % of 
the total emissions. In such a case, the pulsed 
emissions may considerably contribute to the ob- 
served emissions at off-peak orbital phase, and the 
intrinsic ratio of the maximum to minimum fluxes 
of the X-ray emissions from the inter-binary shock 
is larger than F max /F min ~ 2 — 3. This results in 
increase in the flow velocity estimated from equa- 
tion ©. 

2* Although the pulsed X-ray e mission from 
PSRB 1957+20 has been reported bv lGuillemot et al 



1996T ). which corresponds to ~ 0.15 orbital phase. 
Because the emission is concentrated in the for- 
ward direction of the flow, therefore, we expect 
that double peaks due to the Dopplcr effect ap- 
pear at ~ 0.15 phase before and after the phase 
of radio eclipse, which is consistent with the ob- 
servation. 

As we can see in Figure [IJ the observed ra- 
tio of maximum to minimum fluxes is 2 ~ 3. 
For the emissions from the pulsar wind, the out- 
going flux is modified by the Dopplcr effect as 
F„(E) = V 3 F^,(E') and E = VE', with primed 
quantities refereeing to the comoving frame. Here 
the relativistic boosting factor is given by 



(|2012l ). the origin of the pulsed emissions is not 
known. As Table 1 shows, we find that the spec- 
trum average over the whole orbit can be fitted 
better by a power law plus blackbody model. 
The effective temperature and radius of the black 
body component are T e ff ~ 2 x 10 6 K and 
R e ff ~ 0.1 km, respectively, which are explained 
by the core c omponent of t he heated po l ar ca p 



region (c.f. IZavlinl 120071: iTakata et all l2012t) 



V 



1 



17(1-0, cos 0,)' 



(1) 



where r, and /?/ is the Lorentz factor and veloc- 
ity in units of the speed of light of the flow, re- 
spectively, and Of is the angle between the Earth 
viewing angle and the direction of the flow. The 
ratio of maximum to minimum fluxes is estimated 

by 



1-/3/ cosOfj 



F v . m in \1 — Pf COS Of,, 



(2) 



where a is the photon index, 0f, max and Qf.min ~ 
7r — 0f\ max are typical angles at the orbital phases 



The observed flux of the blackbody component 
is Feb ~2x 10" 14 erg cm" 2 s" 1 , which is several 
tens of percent of the total emissions. This flux 
level of the pulsed emission is consistent with the 
result obtained by Guillemot et al. (2012), who 
have estimated the fraction of the pulsed emis- 
sions as about 30 % of the total emissions. 

X-ray tails around pulsars have been inter- 
preted as bow-shocks generated by the supersonic 
motion of pulsars through space, with the wind 
trailing behind as its particles are swept back by 
the pulsar's interaction with the interstellar gas 
it encounters. Studying the diffuse X-ray emis- 
sion may help to better understand their ambi- 
ent environment and the interaction between the 
pulsar and the I SM. Until now only two MSPs, 
PSR B1957 +20 dStappers et aD|2003h and PSR 
J2124-3358 (|Hui fc Beckerl 120061 ). are found to 
be associated with extended X-ray emission. It 
is worth to revisit PSR B1957+20 with a deep 
Chandra observation. A ~ 25" X-ray tail extend- 
ing from PSR B1957+20 with a position angle 
opposite to the pulsar's proper motion direction is 



7 



clearly resolved in the Chandra ACIS image. The 
non-thermal nature and a marginal softening of 
the spectrum of the tail in X-rays as a function of 
the distance from the pulsar support the scenario 
that particles injected at the termination shock is 
dominated by synchrotron cooling. 

Cheng et al suggests that the observed 



length (I) of the X-ray tail can be interpreted as 
the distance traversed by the pulsar within the 
electron synchrotron cooling timescale (t c ), i.e., 
I ~ v p t c , where v p is the proper- motion velocity 
of the pulsar. The cooling time in the X-ray band 



is - 10 8 B m G /2 (Vx/keV) 1/z s, where B mG is the 
inferr ed magnetic field st rength in the emitting re- 
gion ( Cheng et al.l 12006^ . Adopting the inferred 
tail leng th of ~ 9.7 x 10 17 cm at the distance of 
2.5 kpc (|Cordes fe Lazioll2002h and a proper mo 



1/2 



tion velocity of 220 km s 



Arzoumanian et al 



1994 ). the synchrotron cooling timescale is esti- 
mated to be ~ 4.9 x 10 4 yrs. This yields a magnetic 
field of B ~ 17.7 /LtG in the shock region. Con- 
sidering a magn etic field strength of ~ 2 — 6 /xG 
in the ISM (cf. iBeck et ah 20031 and references 
therein), we found that the magnetic field in the 
termination shock might be compressed by a fac- 
tor of ~ 3, which is consistent with th e estimated 
value reported bv lKennel fc Coronitil 
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